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ME/^UREMEHT OF ATMOSPHERIC PRECIPITABLE WATER 


USING A SOLAR RADI(»fETER 

By David E. Pitts » William McAUum, and Alyce E. Dillinger* 
I^mdon B. Johnson Space Center 


SUMMARY 


Measur^ents of total atmospheric water vapor are accomplished hy ratioing 
a water-vaxK>r absorbing region (0.9^3^ micrcmeter) to a clear channel 
(0.8730 micrometer). Aerosol extinction is \isually a slowly varying ft;".ction 
of wavelength in the visible and near infrared. This ratioing effectively elim- 
inates the aerosol and Rayleigh scattering leaving only the water-vapor absorp- 
tion effects. Because a smeill amount of water vapor may satiurate the strong 
water bands and a large amount of water vapor may not sat\u*ate the weak bands, 
this ratio is neither a linear nor an exponential function of precipitable 
water. 


IHTRODUCTIOH 


The purpose of this paper is to describe a method of using a solar radiom- 
eter to deduce precipitable water in the atmosphere. The technique uses two 
channels. One chamiel covers a water-vapor absorption region, and the other, a 
nearby region in which water-vapor absorption is absent. The method operates 
on the €issumption that the aerosol and Rayleigh scattering optical depths are 
approximately the same in the two channels. This assumption allows the ratio 
of the two channels to be used to eliminate the aerosol and Rayleigh scattering 
and thus gives the precipitable water in the slant path. By dividing the water 

1 2 

amount by the air mass, the precipitable water in the verticcd. path may be 
obtained. Such an instznmient was reportedly constructed and calibrated against 
radiosonde data during the early 1960's, but the details of the instrument and 

3 

its success have not been verified. 


•Lockheed Electronics Co., Houston, Texas 77058. 

■^!he air maiSB is defined as the secant of the angle from the zenith. 

^Precipitable water (centimeters) is defined as the total amount of water 
vapor (grams) in a column above a surface area, of 1 squaure centimeter, divided 
by the density of liquid water (ref. 1). 

^Jurcray and Baurker (ref. 2) accredited Foster, Volz, and Poskett with 
the successful construction amd calibration of this instrument. 
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Superscripts : 

- downwelling quantity 

X experimentally determined constant = 0.57^* 

Subscripts ; 
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a atmospheric gaseous absorption processes 
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d\janmy variable 
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RADIOMETER MEASUREMENTS AND ACCURACY 


Atmospheric aerosol and molecular scattering and absorption by atmospheric 
gases are deleterious conditions that distort the upiwelling terrestrial radia- 
tion and cause problems in the automatic classification of multispectral imagery 
by automatic pattern recognition techniques (refs. 3 and k) . Moreover, these 
atmospheric conditions have large temporal euad spatial variations that may re- 
quire atmospheric corrections to be made over a grid covering the field of view 
of the remote sensor. For example, large gradients of precipitable water exist 
across dry fronts in the Midwest. These gradients may cause variations in pre- 
cipitable water of from 1.9 to 3.0 centimeters over a distance of Uo Kilometers 
(ref. l). To allow for such atmospheric variations, a seven-channel Sun radiom- 
eter for measuring the optical depth of the atmosphere was designed, constructed, 
and tested in the field for more than a yea'* (fig. l). This unit is similar to 
the Volz Sun photometer used by Flowers et al. (ref. 5) and Volz (ref. 6) in 
U3 cities in the United States to determine atmospheric turbidity for air qual- 
ity monitoring and is also similar to the multiple wavelength radiometer de- 
scribed by Shaw et al. (ref. 7). To distinguish between different types of 
aerosols (e.g., continental and merit ime) , Reeser (refs. 8 and and Kleen 
(ref. 10) determined that it was necessary to measure atmospheric extinction at 
multiple wavelengths. 

A device was needed to adequately measxire such atmospheric extinction and 
to cover the principal visible and near infrared bands of the Earth Resources 
Technology Satellite ^ERTS) multispectral scanner (MSS) (ref. 11 ) and the Earth 
Resources Experiment Package (EREP) multispectral scamper (ref. 12). An anal- 
ysis (ref. 13) was made that indicated an optimum set of wavelengths 10.38, 

0.50, 0.61, 0.7W79 0 . 8730 , and l.OU micrometers) that could be used to meet 
this need. A study of the accuracy of the operator-radiometer-atmosphere sys- 
tem revealed a ±3 0-percent error in measurements of the aerosol optical depth. 

An additional tone, at 0.9^35 micrometer was chosen to measure absorption caused 
by atmospheric water vapor. 

Kleen (ref. 10) showed that by ratioing the meter reading in the 
0 . 9 ^ 35 -micrometer band to the meter reading in the 0. 8730-micrometer band, the 
water- vapor absorption could be measured. The beuid at 0.8730 micrometer was 
chosen over that at 1.04 micrometers because the former is closer to the peak 
resiKonse of the silicon detector used in the radiometer. A network of these 
radiometers in the field of view of ‘‘he ERTS or EREP multispectral scanners 
will give a small-scale picture of b^ ;h aerosol and water-vapor effects on tra- 
versing signals. However, a more practical operational scheme that uses space- 
craft data for this purpose will have to be developed to facilitate pattern 
recognition of large or remote areas of the world. 


U 



Cousin et al. (ref. Ik) have described the method (the Prediction of 
Response of Earth Pointed Sersors-Reconstruction of Target Reflectance (PREPS- 
ROTAR)) for using these data to correct multi spectral scanner data from ERTS at 
bands of 0.5 to 0.6» 0.6 to 0.7* 0.7 to 0.8, €uid 0.6 to 1.1 micrometers for 
aerosol effects using a table look-^up scheme based on a solution to the radia- 
tive transfer equation known as the doubling method (refs. 15 and 16). With 
the calibration of the water channel of one solar radiometer represented by 
this paper, the other l6 radiometers may be cross-calibrated. The preciritable 
water obtained may then be used as an input, and water- vapor absorption erfects 
may be removed from bands at other wavelengths for ERTS, EREP, or aircraft data 
with a similar look-up table. This look-up table was generated by modifying a 
compressed line-by-line atmospheric transmission model by Deutschman and Calfee 
(ref. 17) by using a set of atmospheric soundings to represent an extreme range 
of atmospheric precipitable water over extremes in temperature and moisture 
distributions (ref. 1). These results will then be included in the PREPS-ROTAR 
correction program. (A paper describing this look-up table for the ERTS MSST 
(0.8 to 1.1 micrometers) and the effect of water vapor on signature extension 
for corn and soybeans is being prepared by the authors as a companion paper.) 


CALIBRATION OF THE SOLAR RADIOMETER 


Atmospheric transmission T(X) for the 0.9^35-micrometer water band was 
calculated by using a least-squares fit of the transmission of the filter fcr 
solar radiometer unit U, as determined from a spectrophotometer (fig. 2). The 
Gaussian function 4>(X) used for the fit was of the fom 


4>(X) = A expiB^X - ^o) J 


Becaxise the spectral response of the silicon detector changes by less than 
0.5 percent and the solar irradiance changes by less thar. 2 percent ever the 
width of the filter function, it was not necessary to include these terms in 

calculating the average transmission f for the beuid irom the atmospheric 
transmission t( A ) . 


T = 


-0.950pm 

I T(X)c>(X)dX 

^.935 


^0.950um 

I (p(X)dX 

•' 0.935 


( 2 ^ 
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The Integrals in equation (2) were evaluated using the Simpson rule, in- 
coi^orating atmospheric transmissions obtained for a 10-layer atmospheric model, 

and calculating water-vapor absorption with about 1 centimeter"^ accuracy 

(fig. 3) smoothed by a triangular filter function 20 centimeters”^ wide at the 
base (fig. U). Valxaes were dalcxilated for 10 radiosonde data sets fref. l) 
representing a range of 0.l66 to 5.851 centimeters precipitable water (fig. 5)> 
This information includes the extreme range of precipitable water observed in 
the United States over a 9-month period in 1972 and 1973. 

The solution for a downwelling known intensity such as the Sun passing 

through an absorbing atmosphere where no scattering occxurs, is 



However, if the atmosphere does not radiate strongly at the wavelengths consid- 
ered and measurements are made at the Earth surface t =0, equation (3) 
simplifies to 


r (0,y) = I e 
va ov 


-T 


o,va. 
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How, if single scattering without an absorptive process is considered, a similar 
expression can be derived as a solution to the radiative transfer equations 
(ref. l8) 
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vs 
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(0,u) = I^^e 
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(5) 


If these two processes are visualized as distinct phenomena with each occurring 
in a separate mathematical layer representative of the real atmosphere where 
both phenomena occur, then the total intensity transversing the layers is 
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Now, if the solar intensity is measxired in two bands, one with aerosol and 
water-vapor absorption and one with only aerosol absorption, if finite filter 
functions and are \ised for each, and if is a slowly varying 

function between to and to Vj^, then equation (5) becomes 
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and equation (T) becomes 
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where to covers the 0.8730-micrometer band and to 

the 0.9^20-micrometer water band. Now, if aerosol optical depth 
slowly varying function of frequency over the range of frequencies 
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-T 

then e can be removed from the integrals in equations (8) and (9) so 

that the ratio of equations (8) and (9) gives 
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where Fg is the filter function over the * Jtter band to v^. Now, since 

ivsa^0’»‘^ = V( 0 . 9 U 20 ) 1^3(0.^) = C2 J(q q^ 3 qj, then equation ( 10 ) becomes 
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where F2(v) » A exp[B(X - Furthermore 
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For wavelerigths near 1.0 micrometer, K is approximately independent of z 
(i.e., not dependent on temperature) so 


X 

o,v,a 


K P P„ n 

V o HgO 


(13) 


8 



so 


C *^(0-9420) 
‘^(0.8730) 



(iU) 


where 0 = — = 
A 

show that it is 


•^0(0.8730) ^ 

*^o(0.9U20) 
of the form 


Numerical evaluations of the preceding integral 
exp(-KP^*m*) so that 


’^0(0.8730) ^(0.9 h 20) 

‘^o(0.9*»20) *^(0.8730) 



(15) 


if the ratio 8730)^'^o(0 9U20) includes the terms C^/A. 

The water-vapor transmission is not an exponential, function of precipita- 
hle water (fig. 5). This fact indicates that the transmission does not obey 
Beer's law. This is because the centers of strong bands saturate at low water 
amounts* leavirtg the weak absorbing lines to contribute to additional absorp- 
tion. The ad hoc function 


T = 



(16) 


was deteimined to represent a good fit to the curve in figure 5 for solar ra- 
diometer unit U. 


Equations (lU) and (13) assume the optical depth fcr aerosol and Rayleigh 
scattering is approximately the same in these adjacent channels. The Rayleigh 
optical depth changes by about a factor of 2 for these wavelengths, but it is 
a very smadl value of approximately 0.01 ± 10 percent, whereas the aerosol op- 
tical depth is about 0.5 (ref. I8). 


Raclosonde data from downtown Houston, Texas, and solar rc.diometer data 
from the NASA Lyndon B. Johnson Space Center (JSC) at Clear Lake City, Texas, 
were collected for the same days. The ratio of J^^ 9l*35)^*^(0 8730) 

pared to the precipitable water amounts from the radiosonde data in fig\ire 5 
and gave a subjective best fit of O.UU for the value of 8730/ 


‘^o(0.9U35)^ 


that represented the curve adequately except at higher water 


Q 



aaK>unts. These anoaalous values are very likely attributable to the higher 
■oisture anounts present at JSC, which is UO kiloaeters southeast of the 
downtoim Houston area where the 0600 l.s.t. radiosonde is launched. Most solar 
radioaeter neasureoents are aade at 0900 to 1100 l.s.t. when the sea breeze is 
starting to penetrate the coastal areas. 

Once the calibration value 8730 )^^ 0(0 9 ^ 35 )^ known for radiom- 
eter unit li, the aeter readings 9 ^ 35 ) "^(O 8730 ) radic.ieter 

unit can be used to obtain the transnissicai in the water band T, which allows 
the precipitahle water in the slant path to be calculated using the curve in 
figure 5* Dividing the precipitahle water in the slant path hy the air mass 
(sec 0 ) will give the precipitahle water in the vertical path. 


Precipitahle wator aeasureaents using othrar radioaeters may be obtained by 
taking the ratio of the measureoents 9 ^ 35)^*^(0 8730 ) nultiplying by 

the calibration value 8730 )^'^o (0 9 ^ 35 ) appropriate radiometer 

froai table 1 . 


Figure 6 will then relate this calculated value to precipitahle water in 
the slant path. The precipitahle *.'ater in the vertical path msf be obtained by 
dividing the slant path amount by the air mass m. 


Because each filter has a slightly different transsdssiou function, the 
resptHise of transmission as compared to precipitahle water is different for the 

0.9^35*^croiBeter band of each radiometer. Figure 6 shows the calibration 
curves for all 17 radiometers. Ihe curves were obtained hy integrating the 
measured filter response functions over the absorption caused hy water atmos- 
pheric vapor; five different atmospheric profiles were used in the same manner 
as Just described for radiometer unit li. Hany of the curves are, for practical 
purposes, identical. 


Once the calibration values 


1-’o(0.8T30)/'’o(0.91.35)\ “* 


radiometer unit k, calibration of the other units that were available were un- 
dertaken by simultaneous measurements with radicmeter unit U and the un- 

known radiometer. The technique used for finding the calibration values 
^^ 0(0 8730)^’*^o(0 9k35)^ nncalibrated radiometers is as follows: 

1. Obtain two sets of simultaneous readings with radicamter unit U and 
the luoknown instrument for the 0.8730- and O.9k35*^crometer channels. The two 
sets should be coonared to eliminate reeding or writing errors. 


2. Use tJo(Q.8730)^'^o(0.9k35)\ 
measured values U(o.9k35)^'^(0.8730)^ 
culate Tj^, where T,^ - ®****^^*^(0.9k35) 


for the calibrated photometer and 

for the calibrated photometer to 
k 

^*^(0.8730) 


the 

cal 
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3. Use figujre 6 to calculate precipitable water in the slant path based 
on oeasurements from radloaeter unit U. 

Use the pz^cipitable water from step 3 to find the transmission for 
the uncalibrated radiometer using figure 6. 

5. Use measured values 9l»35) *^(0 8730) uncalibrated 

radiooteter and the from step U to calculate the calibration constants 

'^o(0.8730)^‘^o(0.9‘»35) * 


^oCO.8730) ^ ^(0»8730) ^ 
*^o(0.9H35) '^(0.9‘»35) “ 


Table I gives calibration values dbtcdned in this manner for all radiameters 
that were available in the susraer of 1973* 

6. Obtain a sec<md set of simultaneous measuraoents on a different day, 
perform steps 1 througb 5» and cross correlate the values of water obtained 
with radiometer unit 4 and the instnmient being calibrated. Figure 7 is an 
example of such a calibration verification. 

If radiometer unit U becomes uncalibrated, units 1, 2, or 3 should be used 
as the standaxd^. 

In severed studies, Morrissey and Brousaides (ref. 19), Barnes et al. 

(ref. 20), and others discuss temperature-induced errors in humidity, problems 
in repetition, and problems with response of radiosondes that served as the 
andard for measurasent of water vapor for this technique. Brouseddes and 
Morrissey (refs. 21 and 22) discuss an improvement in the radiosonde that elim- 
inates the systematic trend to measure drier daytime profiles than infrared and 
other water-vapor measurements in the Barbados Oceanographic and Meteorological 
Experiment project. Their redesign of the radiosonde that was implsaented by 
the Batlooal Oceanic and Atmospheric Administration in 1972 shows an increase 
of 30 i>ercent relative himidity (e.g. , from 30 to 80 percent relative humidity) 

p p 

for some cases above ^00 x lo r/m ($00 millibars). All the data used in this 
report were from the inproved radiosonde syston. 

According to Brousaides (ref. 23), most radiosondes have laboratory accu- 
racies of ±3 percent or less for relative humidity at teispe *atures above 298 K; 
however, at temperatures below 2$3 K, the error increases to ±6 percent with 
a mean bias of as much as 7 percent. However, at present, no assessment of the 
absolute accuracies tar tenperattore and relative humidity is available to cal- 
culate the accuracy of determining total water from radiosonde data. 

Hater vapor is one of the most highly variable of all atmospheric constit- 
uents; thus, the accuracy of callhsrstlag a solar radiemeter by a balloon-launched 
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sensor requiring tens of ndnutes to measure . >rofile will always be limited. 
It is the authors' opinion that the relative snidity measurement from the 
radiosonde is a'^curate to about ±10 percent for values of relative humidity 
gireater than 20 percent. If absolute radiance could be derived for the two 
channels of the radiometer, the accinracy of the radiosonde wovild no longer be 
the limiting factor in determining precipi table water in the atmosphere. 


CONCLUDIBG REMARKS 


I%e technique of using the solar radiometer allows a self-consistent meas- 
uranent of precipitable water to be made that can be used as an input to a trans- 
mission look-up table for other sx>ectral bands of water-vapor absoiption. As 
such, precipitable water serves se a convenient, physically meamingful, inter- 
mediary quantity. 


lyndon B. Johnson Space Center 

RationcLl Aeronautics and Space Administration 
Houstcm, Texsis, February 26, 197*» 
951 - 16 - 00 - 00-72 
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Figtire U.- Calculated transmission through the atmosphere smoothed by a 

triangular filter function 20 centimeters”^ wide. Albany, N.Y., Jan. 
19T2, radiosonde, 0.202 centimeter precipitable water. 
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Photometer 
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O Houston, Feb. 29, 1973 
□ Houston, June 4, 1973 
A Houston, June 2, 1973 
O Houston, May 29, 1973 
V Houston, Nov. 27, 1972 
<1 Houston, Nov. 28, 1972 
I> Houston, Sept. 18, 1973 
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data for 
theoretical 
transmissions 


1 Key West, Fla., July 17, 1972 

2 Lake Charles, La., Sept. 5, 1972 

3 Brownsville, Tex., Oct. 20, 1972 

4 Boothville, La., Nov. 11, 1972 

5 New York, N.Y., Dec. 10, 1972 

6 Tucson, Ariz., June 30, 1972 

7 Albany, N.Y., Jan. 30, 1973 

8 Caribou, Maine, Jan. 10, 1973 

9 Ely, Nev., Dec. 5, 1973 

10 Rapid City, S. Dak., Jan. 10, 1973 
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Figxire 5.- Precipitable water calibration for solar radiometer unit 1*. 
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O Units 1 and 3 
□ Units 2 and 17 
A Unit 4 

O Units 5, 9, 10, 11, 13, and 16 
V Units 6 and 14 
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Figure 6.- Calibration curves for 17 solar radiometers. 
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Figure ?•“ Calibration verification beween solar >"adio.'net er: 



